Phase coherence, inelastic scattering, 
and interaction corrections in pumping through quantum dots 
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Adiabatic quantum pumping in noninteracting, phase coherent quantum dots is elegantly de- 
scribed by Brouwer's formula. Interactions within the dot, while suppressing phase coherence, make 
Brouwer's formalism inapplicable. In this paper, we discuss the nature of the physical processes 
forcing a description of pumping beyond Brouwer's formula, and develop, using a controlled adia- 
batic expansion, a useful formalism to study the effect of interactions within a generic perturbative 
scheme. The pumped current consists of a first contribution, analogous to Brouwer's formula and 
accounting for the remanent coherence, and of interaction corrections describing inelastic scattering. 
We apply the formalism to study the effect of interaction with a bosonic bath on a resonant level 
pump and discuss the robustness of the quantization of the pumped charge in turnstile cycles. 

PACS numbers: 73.23.-b,72.10.-d,73.63.Kv 
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Introduction - Adiabatic quantum pumping in quan- 
tum dots (QD) is one of the simplest, yet non trivial set- 
tings where the interplay between phase coherence and 
non-equilibrium dynamics can be studied both theoreti- 
cally P, i, i, i, H @, 0, H and experimentally Q. In a 
quantum pump, a slow, cyclic variation of the system pa- 
rameters in time generates a direct current through the 
system For suitable cycles the pumped charge may 
be quantized [J, [13], a fact that might have important 
implications for metrological purposes. 

In noninteracting, phase coherent systems, this phe- 
nomenon is encoded elegantly in Brouwer's formula 0, 
describing geometrically the charge pumped per cycle in 
terms of the instantaneous scattering matrices of the sys- 
tem. Pumping in interacting systems is much less un- 
derstood. In the past years several ad hoc approaches 
have been developed to address specific interactiiig prob- 
lems 0. More recently two interesting studies [1, 0, Q 
aimed at deriving generalized pumping formulas, anal- 
ogous to Brouwer's one, within different approximation 
schemes (the average time approximation [y, 0] , and lin- 
ear response 0). The formulas so obtained, while in 
principle equivalent, are both hardly usable as a starting 
point to address generic pumping problems. In addition, 
the lack of a general understanding of the physics of in- 
teracting quantum pumps makes it difficult to develop 
systematic approximation schemes. 

The purpose of this paper is twofold. First of all, 
we develop a systematic physical understanding of the 
new elements brought in by interaction in the physics 
of quantum pumps. Secondly, by making use of a con- 
trolled gradient expansion [J|, we derive a handy for- 
malism to describe the effect of interactions on pump- 
ing within a generic perturbative scheme. The current 
can be conveniently written as the sum of two terms. 
The first, analogous to Brouwer's formula (but expressed 
in terms of the elastic S-matrices), accounts the rema- 
nent phase coherence. The second term express the cor- 



rections due purely to interaction and contains the new 
physics brought in by it: inelastic scattering We 
demonstrate this by studying the effect of the coupling 
to a generic equilibrium environment on a resonant level 
quantum pump [lo| . Computing the dependence of the 
pumped charge on the bath temperature, we prove the 
robustness against interaction of the quantization of the 
pumped charge in turnstile cycles . 

Inelastic scattering and pumping - We begin with a 
qualitative discussion of the effects of inelastic scattering 
on adiabatic quantum pumping. To illustrate our ideas, 
we consider a concrete example of interacting quantum 
pump: a single level quantum dot (QD) coupled to a 
generic bosonic environment in equilibrium at tempera- 
ture T. The Hamiltonian is 

H{t) = ^ efc 4„ Cka + e{t) d'' d + y^^ujqalag + 

k,a q 

+ H [^"(*) Cfc„ + h.c] +dU^Xg (4 + a^g) . (1) 

k,a q 

Here the c?'s refer to the resonant level, the a^'s de- 
scribe the bosonic modes, while Cka annihilates an elec- 
tron in the left/right leads (a = L,R). The chemical 
potential of the two leads is assumed equal and put to 
zero. The instantaneous strength of the coupling to the 
leads is characterized by the parameter T{t) = Fq = 
2 TT v\Va{t)\'^ , where v is the density of states in the 
leads at the Fermi level and the parameters Vq (as well as 
e) are assumed to vary cyclically in time with frequency 
fi. The environment is specified through its spectral den- 
sity 



(2) 



where C is a constant, lUc a high frequency cutoff, and s 
1 corresponds to a sub/super-ohmic bath. We assume a 
small relaxation rate 7 <C F for the bosonic modes. 
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In the adiabatic limit ^ min[r(t),7] the pumping 
current /^(i) from the lead a to the quantum dot can 
be cast as the sum of a first term, which we attribute to 
elastic processes and a second, purely interacting term 
accounting for inelasticity, (t) = (t) + /™ (t) . The 
first 



E 



duj 
2^ 



(-/') Im 



(3) 



is analogous to Brouwcr's formula [2j, though written 
in terms of the instantaneous elastic S-matrix of the 
QD. Indeed, while here f{oj) is the Fermi function of 
the leads, S{uj,t) is defined as the Wigner transform 
S{uj, t) = J dTe'"^''S{t + t/2, t - r/2) of the time depen- 
dent scattering matrix 

S^,li{t,t')^6^fi5{t-t')-i2^uV:{t)G'{t,t')V0), (4) 

where C is the full QD Green's function [l^. When 
f2 <C r,7, we may express using the Fisher-Lee relation 
Sa,i3{io,t) = 5af3 - i2Tr vV*{t) C {u,t)Vp{t) in terms of 
the instantaneous QD Green's function G^{uj,t), solving 
the time independent problem with frozen parameters 
{e{t),V{t)} Following Ref.[il], the latter is easily 

seen to be the expression of the instantaneous elastic S- 
matrix . 

The second term , in its most general form given in 
Eq. p^ - (PT|) . is associated to the effect of inelastic scat- 
tering on the pumped current. This statement is cor- 
roborated by the solution of the present problem at low 
temperatures (T <C F) , where second order perturbation 
theory in the QD-bath coupling applies. Indeed, assum- 
ing for simplicity real QD-lead couplings we obtain 



2t:^ ^ ' T[t) 



Im [To^^9iTo^J , 



(5) 

where the noninteracting T-matrices To'^'" £^re defined by 
the Fisher-Lee relation as 



27TiyV:{t) G',^''iu;,t)Vp{t), 



(6) 



in terms of the instantaneous advanced/retarded QD 
Green's function Go''' (u;,<) = {uj-e{t)T^T{t))-\ Notice 
that /™ contains explicitly the inelastic scattering rate 
associated to the coupling to the environment 



-pin 



de 



[J(e)-J(-e)]A(^ + e,t){/(c. + e) + 7V(e)}, (7) 



where A ~ —2 Im[Go], and iV(e) is the Bose distribu- 
tion function. As a consequence, this term plays no role 
at T = as a result of the identity r™(w = 0) = 0. 
The limit of weak inelasticity/dephasing considered here 
amounts to F™ < F. 

In order to illustrate the role of inelastic processes in 
a concrete pumping cycle, let us consider the turnstile 
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FIG. 1: The turnstile pumping cycle. 

cycle of FiglD keeping the width of the level constant in 
time T{t) — 2r, the pumping parameters are chosen to be 



the energy level e S [— eo, cq] and ST = 



-F,F]. 



This cycle consists in a periodic opening/closing of the 
coupling to the left /right leads, followed by an inversion 
of the position of the level e when the dot is coupled 
to one lead. In the absence of interaction and at low 
temperatures T F the pumped charge is 



Qo = - 

IT 



Arctg(a:;) 



1 



8ttT' 



(8) 



where x — eo/T and the temperature dependence is due 
to thermal broadening only. Notice that, for a noninter- 
acting resonant level, Qo — 1 for x ^ 1. 

Let us now account for interactions and clarify their ef- 
fect on both the elastic and the inelastic channels by sep- 
arating their contributions to the pumped charge. The 
pumped charge associated to Eq.Q is Q'^' = Qq + SQ'^'' 
consists of the sum of the noninteracting result Eq.®, 
plus a correction originating from the loss of unitarity 
of the elastic S-matriccs. This second contribution, ob- 
tained by expanding Eq.([3]) to second order in the QD- 
bath coupling , reads 



6f2 



4a;(9 



^ + 12Arctg(a:) 



,(9) 



T < 
l)(Li, 



F, and 

i-ey) + 



where we set for simplicity 7 ^ 
a, = /dy(3/8 7r2) Sech(f)2{-F(s + 
Us+ii-e-y)) + (3A')r(l + s)Lii+,(l)}, Li,(z) being 
a polylogarithm and T(z) the Gamma function. Notice 
that summing up both contributions we obtain 



CT 



s+l 



F2 



for .T 3> 1. 



(10) 



Intuitively, the growth of the inelastic component 
should be at the expense of the elastic one. In other 
words, SQ'^'' and Q*" should have opposite sign. The 
presence of inelastic channels, not accounted for in the 
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elastic S-matrix, is also behind the lack of quantiza- 
tion of Q'^^ for X ^ 1. Nonetheless, a proper ac- 
count of inelastic scattering processes restores quanti- 
zation for the total pumped charge Q = Q*^' -I- Q™, 
where = § dtP" . Indeed, an computation gives 



where Gj'"' are the full Green's functions of the QD 



r2 



12Arctg(a;) 



Therefore, 



summing up all contributions we have 



Stt 



r2 



(11) 



which is approximately one for a; 3> 1 . This statement is 
valid also in other temperature ranges (e.g. T ^ 7,r). 
Physically, as long as the resonance in the dot can still be 
fully loaded and unloaded (x ^ I) charge quantization 
persists [13] , while for intermediate x the inelastic broad- 
ening of the resonance has a net effect on the temperature 
dependence of Q. 

General formalism - The importance of inelastic chan- 
nels and of keeping track of unitarity pertain also to 
the general pumping formula from which Eq. ^ was ob- 
tained. Let us sketch its derivation. Maintaining full 
generality, we consider a generic multilevel, interacting 
quantum dot described by 

H^^ek cl^ Cka + ^ 4 di + 

+ [v»At)dUka+h.c^ +n,nt[d^,dl]. (12) 

k.a,i 

Here Hint does not involve electrons in the leads. More- 
over, in order to be able to take the adiabatic limit, we 
assume the instantaneous state of the interacting system 
to be Fermi-liquid-like, characterized by an exponential 
decay of instantaneous dot Green's functions with a typ- 
ical rate F. 

The average current flowing through the dot may be 
expressed in terms of Keldysh Green's functions. Indeed, 
using 



«^K.,z(i)(4Wcfca(i)) 

k,i 



(13) 



where Na = J^k'^la'^ka, and using standard tech- 
niques [31 it is possible to show that 



/a(t) =«(/(< -tl)®r^a(il,^) + 

-TUt,h)(^!if{t,-t))-zT<At,t). 



(14) 



Here f(t) is the Fourier transform of the Fermi function 
of the leads and the symbol (S> stands for a convolution 
A{t,ti) (g) B{ti,t') = J dtiA{t,ti)B{ti,t'). The fuU T- 
matrices are defined as 

r:5'<(t,i') = 2^^V:(t) G'-^^^<it,t') Ypit') (15) 



levels 15[, and the boldface notation implies summation 



over level indices. Information on the distribution func- 
tion of the dot, and of its deviation from the equilibrium 
f{uj) is contained in G<. 

A compact formula for the pumped current may 
be obtained expressing in terms of the QD self 
energy S< = S< + S<^, where [S<(i, t')]z,j = 
«Ea,fc^".»W/(efc)K:,, (^')exp[-^efc(^-^')] is the self en- 
ergy associated to the coupling to the leads only, and 
Sint is the one containing interaction (as well as tunnel- 
ing) vertices. Using the relation G*^ = G'' (X> S< ® G" 
we may recast the current ()14|) in the form 



i2ttv{YI G'-) ® S< , ® (G'^V^) , 



(16) 



where the time dependence of the r.h.s is analogous in 
structure to Eq. (fT4|) . 

We are now in the position to take the adiabatic limit 
by performing a gradient expansion 0, [l^ . To this end, 
we notice that we may write lait) = J duj/ {2Tr) I{uj,t), 
where I{w,t) is the Wigner transform of the function on 
the r.h.s. of Eq. (|16p . In order to express it in terms of the 
instantaneous quantities we recall that Wigner transform 
of a convolution has the simple expansion (15.] 

[A(SB]^,r = AB + ^{drAd^B-d^AdrB) + .{l7) 
2,1 

In the adiabatic limit, the expansion may be truncated 
to lowest nonvanishing order. Indeed, each derivative 
produces a factor proportional to Sl/F ^ 1. 

It is crucial to note that, in order to expand Eq. pB]) 
consistently one has to account for conservation of prob- 
ability or unitarity. In equilibrium this is guaranteed by 
the optical theorem, which out of equilibrium generalizes 
to 

np{t, - Kf^it, t') = I j2t:^ ® 



+2 7r^.(V;G'-)®[S^, 



— 1 

int int\ 



(G'^V, 



(18) 



as one can easily see using the Dyson equations. In the 
absence of interaction Eg. p^ reduces to the unitarity 
condition for the time dependent S-matrix, Sa^^ ® 

Slp = 6c.,p6{t-t'). 

Expanding now Ea. ([TC| to first order in the gradients 
one obtains two terms: the first, proportional to the 
Fermi function f{uj), vanishes once we take into account 
consistently the optical theorem Eq. (fTS)) expanded to the 
same order in the gradients. The remaining expression, 
which physically originates from the deviation of the dot 
distribution function from the equilibrium form, takes 
the form 1 = 1^'^+ P" were the first term p'' is given by 
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jm 



J™ , where 



Eq.®, and /™ 

Im[at(V:(t)G-^(u;,t))G'^(c^,t)V„(i)] 



(19) 



(20) 



The equations above, together with Eq.(l3]), reduce the 
computation of the pumped current to that of the instan- 
taneous interacting self energies/Green's functions, ^l^t 
and G''''^, and of their first order expansion in gradients 
(S[^'J'^)'. This quantities may be computed using the 
standard tools of many body theory within a perturba- 
tive scheme. For example, the quantity (S[^"'^) may be 
easily obtained by i ) deriving, e.g. diagrammatically, 
an expression of the full time dependent Til^^^{t,t') in 
terms of the Green's functions, ii ) expanding the latter 
in gradients using Eq. p?| and isolating the term con- 
taining a single time derivative. 

To illustrate this in more detail, let us consider the 
calculation Tiint to second order in perturbation theory 
for the Hamiltonian Eq.([T|). Neglecting the contribution 
of tadpole diagrams, which for <C 7 can be incorpo- 
rated from the beginning in the definition of e(t), the self 

energy reduces to E|,(t,t') = * E, A^cf (t, t')^f (i, *') 
, where are the noninteracting Green's functions 
of the dot, and Df{t,t') = —2Tri{N{ujq) exp{^iu!q{t — 
t')) 4- (1 -I- N{u!q)) exp{±iujq{t — t')) are the bosonic 
propagators. In order to expand these self energies 
in gradients, one needs to first express G^{t,t') = 

GlitM) ® ^tiiiM) ® Gg(t2,t'), where sf (ii,t2) = 
±i2TivY,a^a\ti)f{±ti T t2)V*{t2), and then take the 

Wigner transform of E?((t, t'). Using the formula for 
the gradient expansion and substituting in Eq. p^ - (PT|) 
we find, for real V, Eq.®. 

Relation to previous works - Finally, we notice that 
a direct comparison shows that the formula obtained in 
Ref. 6] corresponds to the sum of I'^^ -I-/™. Therefore, the 
additional term J™ corresponds to the "irreducible ver- 
tex corrections" for the emissivity reported in Ref.0, 17 1 
or to the corrections to the average time approximation 
in Ref.[3]- The formalism presented in this letter re- 
quires the extraction of the self energies with respect 
to a quadratic Hamiltonian. This does not constraint 
future applications to weakly interacting problems only: 
many interesting problems with strong correlations (e.g. 
Kondo effects) can be mapped into quadratic models plus 
deviations from them (e.g. by slave bosons, bosoniza- 
tion, etc.). An alterative approach studying pumping in 
a quantum dot in perturbation theory in the tunneling, 
but accounting exactly for the charging interaction, has 



been recently formulated in Ref. . 

Conclusions - In this paper, we derived a useful for- 
mula for the pumped current through interacting quan- 
tum dots within a perturbative scheme and used it to 
demonstrate that inelastic scattering is the main effect 
forcing a description of pumping beyond Brouwer's for- 
mula. Interestingly, even in the presence of inelasticity, 
pumping can be still described geometrically [sl]. This 
statement draws a clear analogy between pumping and 
the dissipative Berry Phase problem [ll] , whose implica- 
tions, in view of recent work connecting Brouwer's for- 
mula to geometric phases [l9j , are worth of further inves- 
tigation. 
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